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The dynamical properties of a three dimensional model glass, the Frustrated Ising Lattice Gas 
(FILG) are studied by Monte Carlo simulations. We present results of compression experiments, 
where the chemical potential is either slowly or abruptly changed, as well as simulations at constant 
density. One time quantities like density and two times ones as correlations, responses and mean 
square displacements are measured, and the departure from equilibrium clearly characterized. The 
aging scenario, particularly in the case of the density autocorrelations, is reminiscent of spin glass 
phenomenology with violations of the fluctuation-dissipation theorem, typical of systems with one 
replica symmetry breaking. The FILG, as a valid on-lattice model of structural glasses, can be 
described with tools developed in spin glass theory and, being a finite dimensional model, can open 
the way for a systematic study of activated processes in glasses. 
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I. INTRODUCTION 



Upon cooling below the melting point, liquids may ei- 
ther crystallize or enter a super cooled regime. In the 
latter case, as the glass transition temperature Tg is ap- 
proached, molecular motion gets slower and slower and 
the viscosity increases enormously. The relaxation time 
increases by several orders of magnitude and for all prac- 
tical purposes the system remains out of equilibrium. Al- 
though mechanically responding as a solid, structural re- 
laxation is still present, slowing down as the system ages. 
The response to a perturbation applied at a particular 
time tyg will persist for very long times (long term mem- 
ory), preventing the system from reaching equilibrium. 
While the system ages, one time quantities asymptot- 
ically tend to their equilibrium values while two times 
quantities depend explicitly both on the observation time 
and on the time when the perturbation was applied: 
time translation invariance (TTI) is broken, which is a 
manifestation of history dependence. Upon cooling the 
system gets trapped on long lived meta stable states 
which depend on the cooling rate, eventually escaping 
as a result of activated processes. This complex glassy 
dynamics is reflected in several characteristic features, 
as non-exponential relaxation, breakdown of fluctuation- 
dissipation relations, etc, and there is a widespread be- 
lieve that the underlying rugged energy landscape plays 
a fundamental role on the dynamics [|l|. Many aspects 
related to the structural glass transition and the nature 
of the glassy phase are still poorly understood ^^ , as a 
consequence of the inherent complexity of the underlying 
physics and the corresponding lack of a simple but non- 
trivial model. Only recently some analytical progress 
has been made with classic models based on Lennard- 
Jones or soft spheres potentials Q which, on the other 
hand, have been extensively studied via molecular dy- 
namics simulations 1^. Meanwhile, in the field of spin 



glasses, disordered magnetic systems which share many 
physical properties with structural glasses, a reasonable 
theoretical understanding of the basic physics has been 
achieved at least at the mean field level |0,R] . Recently it 
has been found that the equations describing dynamical 
correlation functions of a kind of mean field spin glasses 
simplify, above the transition, to the single equation of 
the Mode Coupling Theory of super cooled liquids [||. 
In particular, several results point out that these mod- 
els (e.g. the p-spin model) with one step of replica sym- 
metry breaking and structural glasses are in the same 
universality class |^Jl^, indicating a deeper analogy be- 
tween the physics of spin glasses and structural glasses 
than previously thought. Nevertheless, it is still mat- 
ter of debate to what extent this analogy can be pushed 
forward. Some important differences between these two 
kinds of systems are evident: while a defining feature of 
spin glasses is the quenched disorder and frustration, no 
obvious quench disorder is present in structural glasses. 
Several models with frustration but no disorder have been 
studied recently ||ll[ . Another important difference is the 
fact that the dynamical variables in a spin glass are lo- 
calized in space and do not diffuse as the molecules in 
a super cooled liquid. We would like to have at hand a 
simple, yet non trivial, microscopic model of super cooled 
liquids in which the fruitful techniques developed in spin 
glass theory could be applied and possibly be extended to 
finite dimension. Moreover, Mode Coupling Theory and 
mean field models in general do not describe activated 
processes, which are responsible for some important char- 
acteristics of the glass phase, as e.g. cooling rate depen- 
dence of the asymptotic state. There are some other very 
simple models with kinetic constraints which reproduce 
quite well the glassy phenomenology [|l2|,|3|, but being 
essentially dynamical they are hard to be studied with 
analytic techniques. 

In this paper we describe salient dynamical features 



of a Hamiltonian lattice model of structural glasses, the 
Frustrated Ising Lattice Gas (FILG). The interesting 
properties revealed by the model rest on an interplay be- 
tween two different kinds of degrees of freedom: diffusive 
particles or translational degrees of freedom and internal 
or orientational degrees of freedom. We present results of 
Monte Carlo simulations of one and two times quantities, 
like density relaxation, density-density correlations and 
responses, correlations and responses of internal degrees 
of freedom, mean square displacements of the particles 
and compression experiments. This three dimensional 
model may serve as a numerical laboratory to perform a 
systematic study of the role played by activated processes 
on the long time dynamics of glass formers, a feature that 
is missing from mean field approximations. 

The paper is organized as follows. In section O we re- 
view the essential results both from equilibrium and out 
of equilibrium dynamics for the FILG. In section |l| re- 
sults for slowly varying chemical potentials are presented 
and compared with the case where a sudden quench is 
performed, while in section IV the breakdown of the 
fluctuation-dissipation theorem in this model is analyzed. 
Finally, in section |y|, we present some conclusions and 
discuss some open questions. 



II. THE FRUSTRATED ISING LATTICE GAS 

The FILG @ is defined by the Hamiltonian: 

H = - J ^ {e^jSiSj - l)ninj - fi^^Ui. (1) 

<ij> i 

The dynamical variables n^ = 0, 1 (i = 1 . . . N) are lo- 
cal densities or site occupations. The Si represent inter- 
nal degrees of freedom, e.g. rotational ones. Although 
molecules in glass forming liquids may assume several 
spatial orientations, here we take the simpler case of only 
two possibilities. Si — ±1. The usually complex spatial 
structure of these molecules is in part responsible for the 
geometric constraints, imposed by the neighborhood, on 
their translational and rotational dynamics. This hin- 
drance effect is mimicked by the quenched random vari- 
ables Eij — ±1. When J — > oo, in order to minimize Ti, 
either the spins should satisfy the bond Ey or at least one 
of the sites i or j must be empty. In this limit, the (site) 
Frustrated Percolation ||lj] model is recovered, where no 
frustrated link can be fully occupied, implying that any 
frustrated loop in the lattice will have a hole and then 
p < 1, preventing the system from reaching the close 
packed configuration. 

The equilibrium properties of the 3D model have been 
studied in |15| for T = 1 and J — 10, which corresponds 
in practice to the Frustrated Percolation limit. In the 
low density regime (/x < 0.75) the behavior is liquid like, 
time correlation functions decay exponentially, equilibra- 
tion is quickly achieved and the particles mean squared 
displacement grows linearly with time, a simple diffusion 



scenario since particles hardly feel any constraint in their 
mobility. At /z « 0.75 there is a percolation transition, 
the corresponding density being p « 0.38. This first tran- 
sition manifests dynamically in the onset of two different 
relaxation regimes in the correlation functions, a fast ex- 
ponential relaxation at short times and a slow, stretched 
exponential relaxation at longer times. While for inter- 
mediate times the diffusion is anomalous, for longer times 
it is linear with the diffusion coefficient becoming smaller 
as the density grows with growing chemical potential (or 
equivalently, with lowering temperature). A second tran- 
sition appears for /j,c ~ 5.5, corresponding to a density 
Pc — 0.67. This is a glass transition at which the re- 
laxation times diverge and the diffusion constant goes to 
zero. In the FILG the structural manifestation of this 
transition is the presence of a frozen percolating clus- 
ter. This transition corresponds to the dynamical one in 
mean field p-spin or Potts glasses or the ideal glass tran- 
sition in the mode coupling theory. Scarpetta et al pq ] 
studied an equivalent, non local version of the FILG (the 
Site Frustrated Percolation, SFP) in 2D, finding a be- 
havior similar to the 3D FILG, the main difference being 
the (Arrhenius) dynamical singularity occurring at zero 
temperature. It was also found in [051 that the spin glass 
susceptibility Xsoi associated with the internal degrees of 
freedom of occupied sites, diverges at the same value of p 
where a glass transition takes place. This is a thermody- 
namic spin glass transition associated with the frozen-in 
of the internal degrees of freedom. On the other hand, 
the compressibility associated with the density variables 
does not present critical behavior but only a maximum as 
observed in many glass formers. In view of this behavior 
one could ask whether a suitably defined nonlinear com- 
pressibility should diverge at the glass transition. This 
should indicate the thermodynamic Al character of the 
glass transition in this model. A nonlinear compressibil- 
ity, HsG = A^(<'?^> — <9n>^), can be introduced with 
qn = N^^ ^^ "^l^l where 1 and 2 label two real replicas 
of the system. One could think that in this quantity, as- 
sociated with fluctuations of an Edwards- Anderson like 
order parameter involving only density variables, a di- 
vergence analogous to the one present in the Xsg might 
be detected. This quantity has been measured [0 but 
no evidence of critical behavior has been found near p^'- 
KsG attains a plateau that seems to be roughly sample 
and size independent. The fact that the spin variables 
present a thermodynamic transition while the particles 
do not points once more to the purely dynamical charac- 
ter of the glass transition. Of course, a thermodynamic 
transition may be present at a lower temperature (or a 
higher chemical potential) corresponding to the Kauz- 
mann Tr- and associated with a configurational entropy 
collapse, but this remains an open question in this model 
which deserves further study. 

At a mean field level, equilibrium properties have been 
studied in closely related models [Q. For large p all 
sites are occupied and the behavior corresponds to the 
Sherrington-Kirkpatrick model, with a continuous spin 



glass transition. Upon decreasing ^ one reaches a tri- 
critical point, below which the transition is first order. 
In this region, analogously to other first order transition 
models like the p-spin and the Potts glass p9|, above 
the critical temperature there is a dynamical transition 
where meta stable states first appear W^. We still do 
not know to what extension this scenario remains valid 
in finite dimensions. 

A characteristic feature of the out of equilibrium dy- 
namics of glassy systems is aging |g^. Slow aging dy- 
namics is present in the FILG both in density-density 
two times correlations as well as in correlations of the 
internal degrees of freedom |pi| after performing sudden 
quenches in ^ from a very small value characteristic of the 
liquid phase (low density) to a high value corresponding 
to the glassy phase. Performing a quench in fj, is similar 
to the application of a sudden compression, we raise the 
density of the system keeping the ratio fi/J < 1 in order 
to prevent close packing. Connected two point correla- 
tions can be defined as: 



C(^6, t^ j 



— ^ n,{t)n^{t^) - p{t)p{t^), 



plateau separates two time scales typical of glassy sys- 
tems: a P (fast) relaxation for short times and an a (slow) 
relaxation at longer times, corresponding respectively to 
the fast movements of the particles inside the dynami- 
cal cages and the large scale, cooperative process that 
takes much more time in order to rearrange the cages. 
Moreover, in this very long time regime {t — t^ ^ t^,), 
the system falls out of equilibrium, the correlations decay 
to zero asymptotically and time translational invariance 
(TTI) no longer holds with the corresponding violation 
of the fluctuation-dissipation theorem (FDT), as will be 
shown in the following sections. A good scaling of the 
correlation functions is obtained assuming a time depen- 
dence of the form h{tw)/h{t + t^) with h{x) given by 
h{x) = exp [(1 — a)~^(x/TY^"'\ where a < 1 and r is 
a microscopic time scale pl| . This form is quite general 
as one recovers the cases of a simple i/i^, dependence 
(full aging) when a = 1 and stationary dynamics when 
a — Q pO]. A similar aging scenario has been found 
associated with the relaxation of internal degrees of free- 
dom ill. 



(2) 



where the global density at time t is given by p{t) — 
N~^'Y^^ni{t). We now define the density autocorrela- 
tions as Cn{t.,tyj) = c{t,tu,)/c{0,tw). 
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FIG. 1. Density autocorrelations after a quench to ^ = 10 
at T = 1 and J = 10 for L = 20. The waiting times range 
from 2^ (bottom) to 2^^ (top) and the averages are over 50 
samples. 

In Fig. n the behavior of C„(i, t^) is shown as a func- 
tion of i in a semi log plot after a quench in chemical 
potential to a value p = 10, for waiting times between 2^ 
and 2^^. A typical aging scenario is present signalling the 
slowing down of the dynamics as the waiting time grows. 
For the longest waiting times the correlation presents a 
rather fast relaxation to a plateau in which the system 
evolves in quasi equilibrium: the dynamics is station- 
ary and the fluctuation-dissipation relations hold. The 



III. SLOW AND FAST COMPRESSIONS 



When the system is gently compressed the glass tran- 
sition depends on the rate of compression . For a fixed 
temperature we have done simulations of slow compres- 
sion experiments, where at each Monte Carlo step (MCS) 
/i is increased by a quantity A/z. Fig. g shows the evo- 
lution of the density for several compression rates, rang- 
ing from 10~^ to 10"^. It can be noted that for slower 
compressions, the final density is higher. Up to jjl k, 1 
the density remains small enough and the system is able 
to equilibrate sufficiently fast, regardless of A/i, and all 
curves collapse. As the density increases, the equilibra- 
tion time also increases since the system cooperative re- 
arrangement involves a larger number of particles in a 
smaller region. In analogy with cooling experiments, for 
each compression rate the system gets out of equilib- 
rium at a different value of p, and the asymptotic den- 
sity increases with decreasing rates. The points where 
the curves for each rate depart from the master curve 
define the glass transition point pg. The dependence 
of /ig on the compression rate is characteristic of glass 
formers and is thought to be consequence of activated 
processes: no cooling (compression) rate dependence is 
observed in mean field models where the glass transition 
appears at a fixed temperature (chemical potential) . For 
a fixed /i, by decreasing the compression rate, the den- 
sity p(/i) tends to the equilibrium density Poo(a') as a 
power law Poo(m) — p(m) oc A/x~^, as can be seen in the 
inset of Fig. 0, the parameters depending on the value of 
/i. The extrapolated equilibrium value is slightly higher 
than the one obtained upon fast compressions, as will be 
shown later. 




FIG. 2. Asymptotic density for different compression rates, 
ranging from 10~^ to 10~® (bottom to top). Inset: power law 
fitting for several values of fi. 



In analogy with temperature cycling experiments in 
spin glasses which show the onset of irreversible pro- 
cesses, here, by increasing and decreasing the chemical 
potential at a given rate, hysteresis cycles can be pro- 
duced whose internal area increases with A/i. Moreover, 
the chemical potential range in which the system behaves 
in an irreversible way also increases with Afi. In Fig. 
the hysteresis curves are presented for several cycling 
rates. Notice that the lower left part of the cycle coin- 
cides with the equilibrium curve and for very slow rates 
the cycle tends to disappear and we recover the above 
annealing equilibrium curve. As before, the maximum 
achieved density increases for slower rates. 
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consider an initial state with an empty lattice, corre- 
sponding to /z = — oo. Then the chemical potential is 
suddenly increased to a value above the critical one. In 
Fig. we show how the density increases after a sud- 
den quench to different /i values, ranging from 1 to 10. 
For all /I < lie — 5.5 the density converges more or 
less rapidly to an asymptotic value, which depends on 
/i. For fi > fic the density increases with a law p(t) 
which seems to be independent of the chemical potential 
and that converges to the asymptotic value with a power 
law: Poo — p(t) oc t~". The best value for the parameters 
are a ~ 0.47 and poo — 0.673, slightly lower than the 
value obtained in the slow compressions. In contrast to a 
sudden compression, by gently annealing the system can 
approach equilibrium more efficiently. Nevertheless one 
can expect that the equilibrium density will be reached 
after a quench in a much longer time scale driven by ac- 
tivated processes. In Fig. ^we show the relaxation of the 
density after a sudden compression to p — 10, deep in 
the glass phase. 
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FIG. 4. Increasing of the density in a L = 40 system after a 
sudden quench in the chemical potential (T = 1 and J = 10). 
All the data for p > 6 collapse on the same curve, which can 
be well fitted by a power law (see text). 
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FIG. 3. Hysteresis curves (pvsp) from cycling experiments 
in a 1/ = 6 system {T — 1 and J = 10). Notice the different 
a;-axis scale of fig. (a). The cycling rates are, from (a) do (d), 
1, 10-\ 10-2 and 10"^ 



Instead of slowly compressing the sample, high den- 
sities can also be obtained by a fast compression anal- 
ogous to temperature quenching experiments. Here we 
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FIG. 5. Decay of the specific volume after a quench to 
fj, — 10 for a lattice with L = 6. The solid line is a power 
law fit over the last two decades, v{t) = t;oo + at~°', with 
a ~ 0.47 and v^ = 1/poo ^ 1.48 (poo ^ 0.673). 



If we assume that the density asymptotic value is es- 
sentially determined by the amount of frustration present 
in a sample, then it is possible to calculate the range of 
allowed values for p^c- In a sample with no frustration 
obviously p^^ = 1, while in a fully frustrated sample 
we have that p™" = 5/8 = 0.625. A fully frustrated 
sample has the interactions Eij such that every plaquette 
is frustrated. The value we find for the random model, 
Poc — 0.673, falls in this interval, as it should, and is not 
very far from the one found in the fully frustrated case. 
This would suggest that the FILG defined on a fully frus- 
trated lattice may be an interesting non random model 
with features similar to those presented in this paper. In 
fact we believe that the glassy behavior of the FILG is 
based much more on its frustration rather than on its 
randomness. 

As long as activated processes do not enter the scene 
we would expect a behavior reminiscent of mean field 
models. In the thermodynamic limit, mean field theory 
for the p-spin spin glass predicts that the internal en- 
ergy will relax with a power law to a threshold value Eth 
greater than the equilibrium one [|lO| . The absence of ac- 
tivated processes will prevent the system from escaping 
the threshold states and it will be kept out of equilibrium 
forever. In a true glass the relaxation of the density at 
sufficiently low temperatures should behave this way. In 
the case of the FILG we do not know what the equilib- 
rium density is, so in order to test the previous hypothesis 
we have measured autocorrelations at constant density. 

Considering again the connected two times autocor- 
relations (eq. 0), we present in Fig. H results of simula- 
tions of quenches at fixed densities. The upper figure 
shows results for a density that is slightly below the crit- 
ical one {p = 0.67) for several waiting times tw There 
exists a characteristic waiting time after which all the 
curves collapse on a master one, meaning that aging is 
interrupted and time translation invariance (TTI) is re- 
covered. The relaxation preserves the characteristic two 
steps with equilibrium dynamics for large i^,. On the 
other hand, for a density p = 0.68, slightly above the 
critical one, the system seems to age forever, as can be 
seen in the lower figure. In this case, even for the largest 
waiting times, the quasi-equilibrium regime is followed 
by an out of equilibrium relaxation where TTI no longer 
holds and fluctuation-dissipation relations are violated. 
The behavior in this two cases is qualitatively different 
and is reminiscent of mean field behavior with a threshold 
or critical value for the density. This will be reflected also 
in the form in which the fluctuation-dissipation theorem 
is violated, as will be discussed below. 
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FIG. 6. Decay of the density autocorrelations after a 
quench at a constant density and L = 20. Note that for 
p = 0.67 (top) and moderate waiting times the dynamics be- 
comes stationary and aging stops, quite different from what 
happens for p — 0.68 (bottom). 

In real space, one way of seeing aging effects is looking 
at the diffusion properties of the particles. The averaged 
mean square displacement (MSD) in the time interval 
(t,tw) is defined by 



R {t,tu 



AT / , ^iV'j l^u 



[n{t)-n{t^)Y 



(3) 



where the function ai{t,tw) is 1 if the i-th particle is 
present at both times t and t^ and is otherwise. In this 
way only particles that do not leave the system enter 
the measures. Obviously, if the density is kept constant, 
fli = 1 Vi. For the smallest waiting times, the density is 
low and diffusion is like in a simple fiuid, R^ being linear 
in t (see Fig. 0). Only for large times the density attains a 
large value, cages are formed trapping the particles, and 
the diffusion is arrested, as signaled by the appearance 
of a plateau. For larger tw and t <^ tw, the particles are 
rattling inside the cages formed by their neighbors and 
the diffusion is still normal, R^ being again linear in t. 
The fact that the first linear regime is longer the smaller 
is iu, is simply because for larger t^ the density at t — 
is larger. As soon as R^ is comparable to the cage linear 
size, the diffusion halts and a plateau is attained. To 



escape from the cage (hence from the plateau), a large, 
cooperative rearrangement must occur among neighbor 
cages. This slow phenomenon is reflected in the exten- 
sion of the plateau. Moreover, for larger times we are es- 
sentially measuring the diffusion of the cages, what can 
be seen as the second linear regime present in Fig. R. 
From this figure it is clear that aging effects with charac- 
teristic dynamical regimes are also present in real space 
properties of the system. 
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FIG. 7. Two-times mean square displacement after a 
quench for J = 10, T = 1 and L = 20. From top to bot- 
tom: smaller to higher t^. 



IV. FDT 

When glasses evolve slowly out of equilibrium 
fluctuation-dissipation relations are no longer valid. In- 
stead, in computer simulations and mean field approxi- 
mations a generalized relation is obtained pO|,pO 22 1: 



R{t, tyj) 



X[C{t,t^)]dC{t,t^) 



T 



dt„ 



(4) 



where C{t,tw) is a two times correlation function and 
R(t,t^) with t > t^ is the associated response. T is 
the heat bath temperature and X{C), the "fluctuation- 
dissipation ratio" (FDR), is a function that measures the 
departure from FDT: at equilibrium X — \ and the usual 
FDT is recovered, while in the out of equilibrium regime, 
X < 1. The function X turns out to be a constant in 
mean field models of glasses, signalling that the asymp- 
totic dynamics of the system is dominated by only one 
time scale (apart from the FDT scale). In this time 
scale, corresponding to the aging regime, the slow de- 
grees of freedom are at an effective temperature differ- 
ent from the bath one and that depends on the FDR 
as T(,ff = T/X |2^. Recently the first attempt at an 
experimental measurement of this effective temperature 
has been done in glycerol glass |2J]. 

In a finite system aging will eventually be interrupted 
due to activated processes leading the system to equili- 
bration at the temperature of the bath. Nevertheless the 



true role of activated processes in the long time dynamics 
of glassy systems still remains to be understood. While 
an analytic approach to the problem is very difficult one 
may gain insight from simulations on finite systems. 
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FIG. 8. Integrated response versus density correlations for 
three waiting times (see text). 

We have measured the FDR on the FILG H. The 
main result is shown in Fig. Q, where X is given by mi- 
nus the derivative of the data curves. In this case we 
have measured density correlations and associated com- 
pressibilities (for details see pq]). This result was ob- 
tained in a lattice of linear size L — 30. The result is the 
one expected for a model having only two dominating 
time scales: the FDT one, when C > Qea and another 
one, when the system is relaxing very slowly out of equi- 
librium and C < (7ea, feAbeing the Edwards- Anderson 
(EA) order parameter 0. The constancy of X in the 
aging regime is exactly what is predicted by mean field 
theory. 

What is very remarkable is that this constancy seems 
to be valid not only for the asymptotic curve {t^ — > oo), 
but also for finite (although large) t^ curves. Then we 
can express the FDR as function only of the smaller time, 
tw, and we observe that it changes very slowly with i^. 

We also note that the fact that X seems to depend only 
on the smaller time is in agreement with some recent ana- 
lytical results. In p^] it has been shown that slow glassy 
dynamics can be viewed as an evolution between "quasi 
states" , where the system (almost) equilibrates before re- 
laxing further. When the field is switched on at time i^ 
it will select a particular set of quasi-states, those with 
the "right" magnetization. If this is the case then the 
externally induced small perturbation ||27| will regress in 
the same way a spontaneous fluctuation does and one can 
show Ig^ that the FDR may depend on the smaller time, 
but not on the larger one. And this is exactly what we 
have observed numerically in the FILG model. 

Nevertheless, in contrast with mean field predictions, 
the curves for different temperatures and chemical po- 
tentials in the glass phase (for a fixed large t^,) remain 
parallel, simply shifting upwards with a corresponding 



decrease in the value of the EA parameter |g5|| . This is 
reasonable as one expects that the FDT regime will be 
the only one present a,s T ^ Tg. This behavior is not 
observed in mean field approximation, where the EA pa- 
rameter has a discontinuous jump at Tg and X changes 
continuously for T < Tg. Other studies in long range 
models seem to point out a continuous approach of Te// 
to the bath one |2^. A systematic study of these ef- 
fects is necessary in order to understand the underlying 
mechanisms responsible for them. 

It is not clear whether the system will undergo also 
a static transition with one step of replica symmetry 
breaking, as the off-equilibrium dynamics would suggest. 
Thermalization at T < Tg is really hard to attain and 
so we did not succeed in answering this question. How- 
ever, even in the case that all the observed dynamical 
features would not correspond to an underlying static 
transition and they would eventually disappear as the 
system equilibrates, the results are interesting for at least 
two reasons. First, real glassy systems can be in the off- 
equilibrium regime all along their life and so it is more 
relevant the study of their out of equilibrium behavior 
than the static one. Second, even if what we are ob- 
serving is due to finite time effects (although very large), 
it is very remarkable that it is well described by mean 
field approximations. Maybe is the case that on time 
scales where activated processes are very rare mean field 
theories (that completely neglect activated processes) are 
really a good approximation. 

We have also measured the FDR associated with the 
relaxation of the spin or internal degrees of freedom. One 
can consider the spin-spin autocorrelation 
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FIG. 9. Integrated response versus spin correlations (see 
text). 



The interesting result which shows up here is that, 
qualitatively, the FDR is the same in the three cases. Fur- 
thermore, fits to the data show approximately the same 
asymptotic value of X. While the data for spin-density- 
spin-density and density-density are exactly on the same 
straight curve the data for spin-spin variables depart from 
the FDT line at a smaller value of the Edwards- Anderson 
parameter. This is due to the fraction of spins placed on 
empty sites. 

The fact that we have measured the same effective tem- 
perature Teff perturbing different degrees of freedom is 
remarkable. It is known that when different degrees of 
freedom are coupled they must reach the same tempera- 
ture. However to our knowledge this is the first time that 
such behavior has been measured in the aging regime of 
a finite dimensional system. 



C,(i,i^) = 1^5,(i)5,(t„) 



(5) 



or, only taking spins from occupied sites, the spin- 
density-spin-density one 



^ns 1,^7 ^w } 



-E 



Si{t)n^(t)Si{t^)ni{t^ 



(6) 



It is worth noting that there is strong evidence |T^] 
that the internal degrees of freedom are responsible for 
a thermodynamic spin glass transition at a temperature 
near Tg, the structural glass transition temperature. The 
constant value of the FDR for the spin variables indi- 
cates that the transition should be typical of systems 
with one replica symmetry breaking. If this is confirmed 
we have at hand a finite model, with short range interac- 
tions, which presents a discontinuous transition, like that 
present in mean field models as the p-spin. In fact, at a 
mean field level, a line of discontinuous transitions have 
been found in models closely related to the FILG [[l8[ . 
More work is needed in order to test this possibility in 
the short range model. 



By applying suitable perturbations, like small random 
magnetic fields, one can define in the standard way the 
integrated responses corresponding to the correlations 
defined above. The parametric plot of the integrated 
responses versus correlations is shown in Fig. 0. 



For the sake of comparison, in Fig. nfl results for 
all three correlations studied are shown for t^ = 10*. 
Notice that the correlation involving both occupation 
and spin variables is always smaller than the other two: 
Cns < 'rnin{CmCs)- For small times, Cns is dominated 
by the particles dynamics, while, for large times, by the 
internal degrees of freedom dynamics. 
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FIG. 10. Several possible correlations Cns, Cs and Cn for 
U = 10^ 



V. CONCLUSIONS 



associated with internal degrees of freedom is observed. 
This may be the first reahzation of a glassy phase char- 
acterized by IRSB in a finite dimensional model. The 
presence of a static transition associated with density 
variables at a low "Kauzniann temperature" Tk < Tg 
remains an open problem in this model. It would be 
interesting to study in more detail how the restoration 
of FDT relations is achieved as Tg is approached from 
below in order to test the strength of the time scales sep- 
aration. Also activated processes may show up for bath 
temperatures larger than those explored in these simula- 
tions as long as frustration begins to be lifted by thermal 
activation. 
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We have described the main characteristics of the out 
of equilibrium dynamics of the three dimensional Frus- 
trated Ising Lattice Gas. It shows a very complex and in- 
teresting phenomenology similar to that observed in real 
glass formers. The glass transition seems to have a purely 
dynamical origin as no evidence of singular behavior was 
found neither in the linear nor in nonlinear compress- 
ibilities. Aging dynamics is present in both density and 
internal degrees of freedom and the decay of autocorre- 
lations is similar to that observed in molecular dynamics 
simulations of Lennard-Jones glasses. Cage effects are 
evident both in density relaxation and on diffusion prop- 
erties where a clear characterization of different regimes 
is possible. From slow compression experiments different 
densities are reached depending on the compression rate 
and hysteresis is observed from cycling experiments. Af- 
ter a sudden compression the density slowly approaches 
a limiting value. Aging experiments at fixed densities 
slightly above and below this limiting value show clear 
qualitative differences showing that this density may sep- 
arate two different dynamical regimes. 

The overall scenario is very reminiscent of mean field 
scenarios for the glass transition, except for the pres- 
ence of compression rate dependence which is absent in 
mean field. Interestingly, the form of the fluctuation- 
dissipation theorem violations is exactly that predicted 
in mean field approximations, i.e. a constant fluctuation- 
dissipation ratio . This means that this model presents 
two very separated time scales. The effective tempera- 
ture in the out of equilibrium scale is nearly temperature 
independent in the glassy phase. This raises the question 
whether activated processes play or not an important role 
in its dynamics. The same constant effective tempera- 
ture is observed associated with the internal degrees of 
freedom. At mean field level, a constant effective temper- 
ature has been obtained in systems characterized by one 
step replica symmetry braking (IRSB). Thermodynam- 
ically, a phase transition with a divergent susceptibility 
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